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Abstract 


The flow field design in bipolar plates is very important for improving reactant utilization and liquid water removal in proton exchange membrane 
fuel cells (PEMFCs). A three-dimensional model was used to analyze the effect of the design parameters in the bipolar plates, including the number 
of flow channel bends, number of serpentine flow channels and the flow channel width ratio, on the cell performance of miniature PEMFCs with 
the serpentine flow field. The effect of the liquid water formation on the porosities of the porous layers was also taken into account in the model 
while the complex two-phase flow was neglected. The predictions show that (1) for the single serpentine flow field, the cell performance improves 
as the number of flow channel bends increases; (2) the single serpentine flow field has better performance than the double and triple serpentine 
flow fields; (3) the cell performance only improves slowly as the flow channel width increases. The effects of these design parameters on the cell 
performance were evaluated based on the local oxygen mass flow rates and liquid water distributions in the cells. Analysis of the pressure drops 
showed that for these miniature PEMFCs, the energy losses due to the pressure drops can be neglected because they are far less than the cell output 


power. 
© 2007 Elsevier B.V. All rights reserved. 
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1. Introduction 


The proton exchange membrane fuel cell is considered to be 
one of the most promising alternative clean power generators for 
portable, mobile and stationary applications because of its low to 
zero emissions, low-temperature operation, high power density 
and fast start-up. In past decades, numerous fuel cell models 
have successfully integrated the electrochemical reactions and 
the transport phenomena [1-23]. 

The flow field design in the bipolar plates is one of key param- 
eters of a PEMFC, which serves as both the current collector and 
the reactant distributor. The reactants, as well as the products, 
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are transported to and from the cell through the flow channels. 
An appropriate flow field design can enhance the reactant trans- 
port and the liquid water removal from the cell to minimize the 
concentration polarization. To this end, various flow field config- 
urations, including parallel, serpentine, interdigitated and many 
other combined versions, have been developed. Many efforts 
have been devoted to the optimal flow field design for obtain- 
ing high cell performance [11,25-38]. Nguyen [24] designed 
the interdigitated flow field which forces the reactants to flow 
through the gas diffusion layer (GDL) with the shear force of this 
gas flow helping to blow out the liquid water that is entrapped in 
the inner layers of the electrodes. As a result, the mass transport 
of the reactants from the flow channel to the inner catalyst layer 
can be improved and the water flooding at the cathode can be 
significantly reduced. Cha et al. [25] investigated the PEMFC 
flow channel scaling behavior with interdigitated, serpentine and 
spiral interdigitated flow fields. Their results indicated that for 
most flow pattern archetypes, the optimal feature size occurs at 
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Nomenclature 


chemical activity of water vapor 

cathode cross-sectional inlet flow area (m7) 
reaction area m?) 

reference exchange current density (A m7?) 
source term of variable & 

mass fraction 

quadratic drag factor 

hydraulic diameter of the flow channel (m) 
mass diffusivity (m? s7!) 

Faraday constant (96,487 C mol`!) 

current density (A m7?) 

current density (A m73) 

coefficient of water vapor condensation rate (s7 1) 
coefficient of water vapor evaporation rate (s~!) 
permeability (m7) 

distance along the channel measured form the 
inlet (m) 

molecular weight (kg mol~!) 

pressure (atm) 

mass flow rate (kg s7! m7?) 

universal gas constant (8.314 J mol`! K7!) 
volume ratio occupied by liquid water 

time (s) 

temperature (K) 

x direction velocity (m s7!) 

y direction velocity (m s7!) 

fuel velocity at the cathode inlet (m s7!) 

z direction velocity (m s7!) 

cathode pressure drop loss (W m7?) 

width of rib (m) 

x direction coordinate (m) 

y direction coordinate (m) 

z direction coordinates (m) 

species valence 


Greek symbols 


electrical transfer coefficient in forward reaction 
electrical transfer coefficient in backward reaction 
porosity 

dependent variables 

phase potential function (V) 

over-potential (V) 

water content in membrane 

viscosity (Pa s) 

kinematic viscosity (m2 s7!) 

density (kg m~?) 

electric conductivity of membrane 

tortuosity of the pores in the porous medium 
exchange coefficient 


Superscripts 


reference value 


Subscripts 

a anode 
aver average 
c cathode 
eff effective 


Ht hydrogen ion 
H2 hydrogen 
H20 water 


in inlet 

k kth species of the mixture 
O2 oxygen 

sat saturation 


total total 


an intermediate channel dimension. Liu et al. [26] investigated 
the use of baffle-blocked channels to improve reactant transport 
and cell performance of PEMFCs. Their results indicated that 
the blockage effects of the baffles forced more reactants through 
the GDL to enhance the chemical reactions to then augment the 
cell performance. Shimpalee et al. [27] investigated the effect of 
the number of gas paths on the reactants transport and cell perfor- 
mance of a 200 cm? PEMFC with serpentine flow fields. Kumar 
and Reddy [28] investigated the effect of the flow field design on 
the steady and transient state performance of PEMFCs with ser- 
pentine, parallel, multiple parallel and discontinuous flow fields. 
Yan et al. [29] investigated the effects of the flow-field design 
(parallel flow field, Z-type flow field, serpentine flow field, par- 
allel flow field with baffles and Z-type flow field with baffles) 
and air flow rate on the cell performance and pressure drop. Yan 
and co-workers [30,31] proposed a novel straight flow channel 
with tapered height and width to improve the efficiency of reac- 
tant utilization in the PEMFCs. Their results demonstrated that 
the tapered flow channel increases the reactant velocity which 
enhances the reactant transport through the porous layers, the 
reactant utilization, and the liquid water removal. The results 
also show that the cell performance can be improved by either 
decreasing the height taper ratio or increasing the width taper 
ratio. Shimpalee and Van Zee [32] investigated the effects of the 
channel and rib cross-sectional areas on the reactants distribution 
and cell performance of 25 cm? PEMFCs with serpentine flow 
fields. They showed that narrower channels with wider rib spac- 
ings give higher performance. Wang et al. [33,34] investigated 
the effect of the flow channel area ratio on the cell performance 
and local transport characteristics for PEMFCs with parallel and 
interdigitated flow fields. Wang et al. [35] proposed a novel 
serpentine-baffle flow field design to improve the cell perfor- 
mance compared that for the conventional serpentine flow field 
design. They used a three-dimensional model to analyze the 
corner effect and under-rib convection effect in the cell and 
found the both effects enhance the oxygen utilization and liq- 
uid water removal capability. Yang and Zhao [36] investigated 
experimentally the effect of the anode flow field design on the 
performance of direct methanol fuel cells, their results show that 
both open ratio and flow channel length had significant influ- 
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ences on the cell performance and pressure drop. An open ratio 
of around 50% led to the best cell performance at moderate and 
high methanol solution flow rates. However, at low methanol 
solution flow rates, a larger open ratio yielded a higher power 
density at higher current densities, corresponding to the mass 
transfer limitation region. In addition, it was found that a longer 
flow channel usually gave better cell performance, but caused a 
larger pressure drop. Xu and Zhao [37] proposed a new flow 
field design termed the convection-enhanced serpentine flow 
field (CESFF), for electrolyte-based fuel cells, obtained by re- 
patterning the conventional serpentine flow fields (SFFs). Their 
experimental results showed that the new flow field resulted in 
substantial improvements in both cell performance and operat- 
ing stability as opposed to the conventional serpentine flow field 
design. 

The design parameters in the bipolar plates for the SFF 
design mainly include the active area, cross-sectional shape 
of the flow channels, widths of the flow channels and ribs, 
heights of the flow channels and ribs, number of serpentine 
flow channels and number of flow channel bends. This liter- 
ature reviews indicates that although there have been some 
analyses of SFF designs, a systematic, comprehensive study 
of the effects of the SFF design parameters is still needed. 
For example, to reduce the computational effort, some stud- 
ies only simulated a symmetrical unit of the cell which might 
be reasonable for parallel flow fields, but is not acceptable 
for the SFF design because the reactant flow in each chan- 
nel differs and the reactants experience different pressure drops 
and concentration changes along the channel. For example, for 
the single SFF shown in Fig. 1, the channel cross-section is 
small, while the channel length is very long. Then, the pres- 
sure drop at neighboring locations between adjacent channels, 
for example, points a and b in Fig. la, would be substantial, 
thus there is a significant pressure gradient across the porous 
electrode, much larger than the pressure gradient along the 
channel direction, resulting in considerable cross-leakage flow 
between adjacent channels (under-rib convection) [35,37,38]. 
Therefore, full-cell models are needed to simulate the flow 
and reactions in the SFF design. This paper describes a study 
of the flow field design parameters for miniature PEMFCs 
with the SFF using a three-dimensional, full-cell numerical 
model. The study focuses on the effects of the number of 
flow channel bends, number of serpentine loops and the flow 
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channel width ratio on the cell performance, with assump- 
tions of a constant active area and flow channel/rib height 
and a rectangular flow channel cross-section. The local oxy- 
gen mass flow rates, liquid water distributions and local current 
densities in the cells for various flow fields are examined. 
The pressure drops are also evaluated for the various SFF 
designs. 


2. Flow channel designs 


This paper focuses on miniature PEMFCs. In general, water 
flooding occurs on the cathode electrode of the PEMFCs, 
because the electrochemical reaction on the cathode produces 
water vapor. If the partial pressure of the water vapor is higher 
than the saturation pressure, the water vapor condenses to form 
liquid water. When a large amount of liquid water accumu- 
lates in the porous layer pores, the oxygen transport resistances 
increase and oxygen mass flow rates decrease. Therefore, the 
cathode flow field design is a key factor for enhancing reactant 
and product transport and removing liquid water. Therefore, 
the cathode flow field designs are studied in this paper. The 
anodes for all cells use parallel flow field designs, with 12 flow 
channels and 11 ribs all with the same heights of 1mm and 
the same widths of 1 mm, while the cathodes use serpentine 
flow field designs with the various design parameters. All cells 
have rectangular channel cross-sections with 23 mm x 23 mm 
active areas and the same thicknesses for the GDL, catalyst 
layer (CL) and membrane (MEM) geometries. All cells have 
the same flow channel and rib heights of 1 mm. The operating 
conditions must remain constant for a fair comparison for all 
cells. The cell temperature is assumed to be 323 K, the reactants 
on the anode side include hydrogen and water vapor with a rela- 
tive humidity of 100%, the reactants on the cathode side contain 
oxygen, nitrogen and water vapor with a relative humidity of 
100%, the inlet flow rate on the anode side is 260 cm? min7!, 
the inlet flow rate on the cathode side is 700cm? min™! and 
the inlet pressures on the anode and cathode sides are both 
l atm. 

The effect of the number of flow channel bends was con- 
sidered using cell designs with single SFFs with 7, 11 or 15 
bends as shown in Fig. 1. To eliminate the effect of the flow 
channel width ratio, the channel and rib widths were equal for 
each cell, e.g. channel and rib widths were 1.533 mm for the cell 


(a) 


(b) 


(c) 


Fig. 1. Schematics of single loop flow fields with various numbers of flow channel bends: (a) 7 bends; (b) 11 bends; (c) 15 bends. 
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A 
(a) 


(b) 


(c) 


Fig. 2. Schematics of single, double and triple SFFs. (a) Single serpentine; (b) double serpentine; (c) triple serpentine. 


Table 1 

Physical dimensions of the fuel cell models 

Quantity Value 
Reaction area (mm?) 23 x 23 
Channel height (mm) 1 

Rib height (mm) 1 
Anode GDL thickness (mm) 0.4 
Anode CL thickness (mm) 0.005 
Membrane thickness (mm) 0.035 
Cathode GDL thickness (mm) 0.005 
Cathode CL thickness (mm) 0.4 


with 7 bends, 1 mm for the cell with 11 bends and 0.742 mm 
for the cell with 15 bends. The effect of the number of ser- 
pentine loops was analyzed for single, double and triple SFF 
designs. The channel and rib widths were | mm for all three 
cells, which included 12 channels and 11 ribs as shown in Fig. 2. 
The effect of the flow channel width ratio was considered for 
cells with the triple SFF design. The flow channel width ratio, 
n, was defined as the ratio of the total flow channel width to the 
total cell width with values of 0.4, 0.522 and 0.6, as shown in 
Fig. 3, where 7 =0.522 denotes that the flow channel and rib 
widths are equal (1 mm). Besides these design parameters, all 
the other parameters, including the thickness, porosity, etc., of 
the GDL, CL and membrane were the same for all cells. The 
detailed physical dimensions of the fuel cells are summarized in 
Table 1. 


(a) 


(b) 


3. Numerical model 


A three-dimensional model was used for the full fuel cell 
to analyze the electrochemical reactions and transport of the 
reactants and products. The cell was divided into the anode flow 
channels, membrane electrode assembly (MEA, including the 
anode GDL, anode CL, proton exchange membrane, cathode CL 
and cathode GDL) and the cathode flow channels. The governing 
equations include the mass, momentum, species and electrical 
potential conservation equations. The model assumes that the 
system is three-dimensional and steady, the reactants are ideal 
gases, the system is isothermal, the flow is laminar, the fluid is 
incompressible, the thermophysical properties are constant, and 
the porous layers such as the GDL, CL and proton exchange 
membrane (PEM) are isotropic. 

The transport equations for the three-dimensional PEMFCs 
are given as: 

Continuum equation 


ðu av 


ðw 


gu B Peg Sc 1 
əx ody az (1) 
Momentum equation 
E F ge T 
u— —+w— 
A ax ay ag 
Eeff OP Pu u Pu 
=s S 2 
p ox + VEeff (Fa ay2 3z2 + Su (2) 


(c) 


Fig. 3. Schematics of triple SFFs with various flow channel width ratios. (a) n =0.4; (b) n =0.522; (c) n=0.6. 


X.-D. Wang et al. / Journal of Power Sources 175 (2008) 397—407 


Eeft OP v v av 
= -AEE y ver (Sat at aa) tS (3) 


fet OP ew w Pw 
= —-— — S 4 
sa een (53 ae a + Sw (4) 
Species equation 
aCk aC; OCx 
Eeff (« Ox +v Jy +w a ) 
PCE PCE 3Cr 
=D Sot S 5 
pat ( zZ iy? + aa ) + Se + SL (5) 


In the momentum equation, £ is the porosity and Su, Sy and Sw 
are the corrected terms of the reactant flow in the gas diffusion 
layer and the catalyst layer and of the proton transfer in the 
PEM, which is listed in Table 2. In the species equation, Dj cfr 
is the effective diffusion coefficient and Se represents the source 
terms due to the chemical reaction in the catalyst layer and the 
proton exchange membrane. The Bruggeman correction [23] 
is employed to describe the influence of the porosity on the 
diffusion coefficient 


Dy ett = Dye* (6) 


If the partial pressure of water vapor is greater than the satura- 
tion pressure, the liquid water forms in the PEMFC. To consider 
the effect of the liquid water formation, itis assumed that the pore 
in the porous material is blocked by liquid water, which results in 
the modification of the diffusion coefficient and the porosity in 
the species equation. The source term Sy, in the species equation 
due to the liquid water is determined by [14] 


Eeff CHO 
Mpy oke. ———— 
SL — H20%c pRT 


Ke€ettS( Psat — Pp,o), 


(Pm o — Pat), if Pmo > Prat (7) 


if Po < Peat 


where the saturation s is the ratio of the liquid water volume to 
pore volume in the porous material, M the molecular weight of 
water, ke the condensation rate constant of water, ke the evap- 
oration rate constant of water and £eff is the modified effective 
porosity of the porous medium by considering the liquid water 
effect which is given by 


Eeft = e(1 — 5) (8) 
The saturation pressure of water can be expressed as 


2 = -572 -Ira 
Poat =10 2.1794+0.02953T—9.1837 x 107? T4+1.4454x 107‘ T (9) 


Table 2 


Source terms in the governing equations 
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To calculate the local current density, the phase potential 
equations should be solved 


0 oP 0 op ð oP 
= S; 10 
we (ae) + ay (a5) ta (on) pi 


where S;=0 in the membrane, $j=—ja in the anode catalyst 
layer, S; =j in the cathode catalyst layer, ® the phase potential 
and om is the ionic conductivity of the membrane. In this study, 
Butler—Volmer equation [7] is used to calculate the transfer cur- 
rent density generated by the electrochemical reaction 


Cy 1 
i — A jref ~H e(%aF/RT)n = 1l 
Ja = AJo,a (%) | ele F/RT)n (1) 
-C 1 
> 4 cref O2 (@aF/RT)n _ 
je = Aji ( c) k a n sans (12) 
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where A j is the reference exchange current density, œa and a 


the transfer coefficient in anode and cathode catalyst layers, n 
the over-potential, R the gas constant and T is the temperature 
of the fuel cell. The ionic conductivity of the membrane can be 
calculated by [2] 


Om(T) = ot! exp | 1268 ee (13) 
= i 303 T 

and the reference ionic conductivity is 

ott = 0,005139A — 0.00326 (14) 
__ J 0.043 + 17.81a — 39.854? + 36.0, O<a<1 a3) 
~ | 14+1.4a-— 1), l<a<3 


where a is the activity of water vapor at the cathode side. The 
relationship between the phase potential and the current density 
is given by 


jpeg (16) 
Ox 
op 
Dai (17) 
e oy 
op 
iz = —Om— (18) 
Oz 


Boundary conditions at the anode flow channels and the cath- 
ode flow channels are as follows: the inlet flow rates are constant, 
the inlet gas compositions are constant and the flows are fully 
developed at the outlets of the anode and cathode flow channels. 
At the solid walls, no slip and zero fluxes are hold. At the inter- 
faces between the gas channels, the diffuser layers, the catalyst 
layers and the PEM, equalities of the velocity, mass fraction, 
momentum flux and mass flux are applied. More details were 
given elsewhere [31]. Main parameters used in the model are 
listed in Table 3. 

The model uses non-uniformly distributed elements with 116, 
93 and 33 elements in the x, y and z directions (the flow channel 
inlets are along the y direction, and the cell height is along the z 
direction). The grid independence was examined in preliminary 
test runs. Four non-uniformly distributed grid configurations 


Table 3 

Fuel cell parameters 

Parameter Value 

Ajg! (Am) 9 x 108 
Ajgt (Am™) 1.5 x 10° 
Qa 0.5 

Me 1.5 

€channel 1 

Tchannel 1 

Kchannel (m?) o0 

EGDL 0.5 

TGDL 1.5 

kapl (m?) 1.76 x 10718 
ECL 0.4 

TCL 1.5 

ker (m*) 1.76 x 107! 
EMem 0.28 

TMem Dagan model 
kMem (m?) 1.8 x 10718 


were evaluated for PEMFC with the parallel flow channel design. 
The numbers of elements in the x, y and z directions were: (I) 
70 x 70 x 25, (II) 93 x 93 x 33, (III) 116 x 93 x 33 and (IV) 
116 x 116 x 41. The influence of the number of elements on 
the local current density distributions at the operating voltage of 
0.3 V for the fuel cell with parallel flow field is shown in Fig. 4. 
The local current density distribution for grid (I) differs from 
those for grids (II-IV) with the deviation about 3.9%. However, 
the local current density distributions for grids (II), (II) and (IV) 
do not show any notable differences. The difference between the 
local current densities for (II) and (IID is about 0.32% and the 
difference between (III) and (IV) is about 0.25%. Grid (IIT) was 
chosen for the simulations as a tradeoff between accuracy and 
execution time. 

In the present paper, the oxygen mass flow rates and liquid 
water distributions on the cathode GDL-CL interface at the two 
representative operating voltages of 0.3 and 0.7 V are analyzed 
to understand the effects of the various flow field designs on the 
cell performance. For cells with various numbers of flow channel 
bends and cells with various numbers of serpentine loops, the 
local transport phenomena in the entire cell were analyzed by 


—o— 70x70x25 
—o— 93x93x33 
—4— 116x93x33 
—v— 116x116x41 


15000 


0.000 0.005 0.010 


y (m) 


0.015 


Fig. 4. Influence of number of elements on the local current densities on middle 
cross-section in the PEM at the operating voltage of 0.3 V. 
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plotting all local physical quantities along a specified serpentine 
channel from the inlet to the outlet. For the double SFF the data 
was plotted along the AA’ serpentine loop while for the triple 
SFF, the data was plotted along the BB’ serpentine loop (as 
shown in Fig. 2b and c). However, when analyzing the effect of 
the flow channel width ratio, the data was plotted at y= 11.5 mm. 


4. Results and discussions 
4.1. Effect of number of flow channel bends 


Fig. 5 presents the /—-Voey polarization curves and I—Weel 
power density curves for single SFFs with the various number 
of flow channel bends. When the operating voltage is greater 
than 0.7 V, the /-V.., curves and J—W-e curves for all cells with 
the various numbers of flow channel bends almost coincide, indi- 
cating that the cell performance is not dependent on the number 
of bends because at the higher operating voltages the electro- 
chemical reaction rate is lower with a small amount of oxygen 
consumption and liquid water, thus the cells with the various 
numbers of bends all provide sufficient oxygen to satisfy the 
needs of the electrochemical reaction. As the operating voltage 
decreases, the current density gradually increases and the cell 
performance for the various numbers of bends starts to differ. 
As the number of bends increases, the limiting current density 
and cell output power density increase and the cell performance 
improves. 

Fig. 6 presents the local oxygen mass flow rates across the 
cathode GDL-CL interface at the operating voltages of 0.3 and 
0.7 V for the various number of bends, where L denotes the dis- 
tance measured from the flow channel inlet. Fig. 6 shows that at 
the high operating voltage of 0.7 V, the oxygen mass flow rates 
are almost unchanged along the flow direction and are almost 
equal for the three different numbers of flow channel bends. 
Since the oxygen mass flow rate on the cathode GDL-CL inter- 
face is equal to the oxygen consumed by the electrochemical 
reaction in the CL, then the number of flow channel bends has 
little effect on the cell performance at high operating voltages. 
At the low operating voltage of 0.3 V, the oxygen mass flow rates 
are all greater than those at the high operating voltage of 0.7 V for 
the three different numbers of flow channel bends, indicating that 
at the low operating voltages the electrochemical reactions are 


--O-- 7 bends 
--O-- 11 bends 
--4-- 15 bends 


Ven (V) 


—o— 7 bends 
—o— 11 bends 
—4— 15 bends 
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Fig. 5. Effect of number of flow channel bends on cell performance. 
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Fig. 6. Oxygen mass flow rates on cathode GDL-CL interface at operating 
voltages of 0.3 and 0.7 V for single loop flow fields with 7, 11 and 15 bends. 


stronger with more oxygen consumption. For all three numbers 
of flow channel bends, the maximum oxygen mass flow rates 
occur at the inlet and gradually decrease towards the outlet, with 
local peaks at the bends of the flow channels. The inlet oxygen 
concentration is the highest, thus the largest oxygen mass flow 
rates through the GDL-CL interface will occur near the inlet. 
Then, the oxygen is gradually consumed by the electrochem- 
ical reaction, so the oxygen concentration gradually decreases 
along the flow direction and, thus, the oxygen mass flow rate 
across the GDL-CL interface gradually decreases. At the chan- 
nel bends, the corners of flow channels force the reactants to 
flow into the GDL and CL; thus, the oxygen mass flow rate 
through the GDL-CL interface sharply increases at the bends. 
Fig. 6 also indicates that at the same L, the oxygen mass flow 
rate increases as the number of flow channel bends increases 
because cells with larger numbers of bends have smaller flow 
channel cross-sectional areas; thus, at the same reactants inlet 
mass flow rates, the reactants inlet velocity higher. Furthermore, 
a larger number of flow channel bends means that there are more 
corners in the channels, which forces more oxygen flow into 
the GDL to increase the local oxygen mass flow rate. There- 
fore, increasing the number of flow channel bends enhances the 
oxygen utilization efficiency. 

The electrochemical reaction in PEM fuel cells produces 
water vapor, which normally condenses since the normal oper- 
ating temperatures are between 40 and 80°C. However, liquid 
water then increases the resistance to fuel transport into the cata- 
lyst layer, resulting in reduced cell performance. Fig. 7 presents 
the liquid water distributions on the cathode GDL-CL interface 
at operating voltages of 0.3 and 0.7 V for the various numbers 
of flow channel bends. The results indicate that the liquid water 
concentration at the high operating voltage of 0.7 V is less than 
at the low operating voltage of 0.3 V because at the high oper- 
ating voltage of 0.7 V the electrochemical reaction is weaker 
with less production of liquid water. At the low operating volt- 
age of 0.3 V, the liquid water concentrations gradually increase 
along the flow direction, while reach the local minimum values 
at the bends of the flow channel. Although more liquid water 
is produced in the porous layers near the flow channel inlets 
due to the stronger electrochemical reaction rates, the high local 
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—o—7 bends, V=0.3 V 
—o— 11 bends, V=0.3 V 
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—e*— 11 bends, 
—«— 15 bends, 


Fig. 7. Liquid water concentrations on cathode GDL-CL interface at operating 
voltages of 0.3 and 0.7 V for single loop flow fields with 7, 11 and 15 bends. 


oxygen mass flow rates in the porous layers generate high shear 
forces which carry the liquid water downstream to the flow chan- 
nel, so the liquid water concentrations gradually increase along 
the flow direction. However, at the bends of the flow channels, 
the corners force the reactants to enter the porous layers which 
enhances the liquid water removal; thus, local minimums appear 
in the liquid water distributions at the bends. At low operating 
voltages, the electrochemical reaction produces large amounts 
of liquid water. If the liquid water is not removed but accu- 
mulates in the pores of the GDL and CL, the oxygen transport 
resistances increase, which reduces the electrochemical reaction 
rate. Fig. 7 shows that the liquid water concentrations with 15 
flow channel bends is much lower than for 11 and 7 flow chan- 
nel bends, because it has a higher oxygen mass flow rate. Thus, 
a proper flow field design can effectively enhance the reactant 
transport rates and liquid water removal, which both improve cell 
performance. 


4.2. Effect of number of serpentine loops 


Fig. 8 shows the J—V¢eq, polarization curves and J—Wee| power 
density curves for the single, double and triple loop flow fields. 
Fig. 8 shows that as for operating voltages greater than 0.7 V, 
the number of serpentine loops has little effect on the cell per- 
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Fig. 8. Effect of number of flow channel serpentine loops on cell performance. 


formance. As the operating voltage decreases, the number of 
serpentine loops begins to affect the cell performance with the 
single SFF design having better performance than the double 
and triple flow field designs. 

Fig. 9 shows the oxygen mass flow rates on the cathode 
GDL-CL interface at the operating voltages of 0.3 and 0.7 V 
for the various numbers of serpentine loops. Fig. 9 shows that 
at the higher operating voltage of 0.7 V, the oxygen mass flow 
rates are almost equal and unchanged along the flow direction 
for the single, double, and triple SFF designs, so the cell perfor- 
mance is not dependent on the number of serpentine loops. At 
the lower operating voltage of 0.3 V, the maximum oxygen mass 
flow rates for all three flow field designs all occur at the inlet 
and then gradually decrease towards the outlet. The oxygen mass 
flow rates have local maximums at the bends in the flow channels 
for the same reason as given in Section 4.1. Fig. 9 also shows 
that at the same L, the oxygen mass flow rates increase with 
decreasing number of serpentine loops because as the number of 
loops decreases the reactant inlet velocity increases. Compared 
with the parallel flow field configuration where the reactants 
only enter the porous layers by diffusion, in the loop design the 
reactants are also forced into the GDL by the flow field design 
which creates under-rib convection. The under-rib convection 
can significantly enhance the local oxygen mass flow rates in 
the porous layers under the ribs. For the single SFF design, 
the under-rib convection occurs between all the ribs, while for 
the double and triple SFF designs, under-rib convection only 
occurs at the specified ribs. For example, under-rib convection 
occurs at rib | in Fig. 2b, but not at rib 2 in Fig. 2b. There- 
fore, the reactant utilization for the single serpentine flow field 
is higher than for double and triple serpentine flow field 
designs. 

Fig. 10 shows the liquid water distributions on the cathode 
GDL-CL interface at the operating voltages of 0.3 and 0.7 V for 
the various numbers of serpentine loops. The liquid water distri- 
butions are opposite to the oxygen mass flow rate distributions, 
as discussed in Section 4.1. The single flow loop enhances the 
reactant utilization and liquid water removal, so it has the best 
cell performance. 
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Fig. 9. Oxygen mass flow rates on cathode GDL-CL interface at operating 
voltages of 0.3 and 0.7 V for single, double and triple loop flow fields. 
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Fig. 10. Liquid water concentrations on cathode GDL-CL interface at operating 
voltages of 0.3 and 0.7 V for single, double and triple loop flow fields. 


4.3. Effect of flow channel width ratio 


Fig. 11 shows the -Vee polarization curves and I—Weel 
power density curves for the various flow channel width ratios. 
The results show that for operating voltages greater than 0.7 V, 
the effect of the flow channel width ratio on the cell perfor- 
mance can be neglected. At the lower operating voltages, the cell 
performance improves slightly as the flow channel width ratio 
increases, indicating that the cell performance is not strongly 
related to the flow channel width ratio for PEMFCs with the 
SFF design. For the parallel flow field design, as the flow channel 
width ratio increases the cell performance significantly improves 
because the reactants are transported into the GDL and CL 
mainly by diffusion. A larger flow channel width ratio increases 
the contact area between the reactants and the GDL, which 
allows more reactants to directly diffuse into the porous layers 
to participate in the electrochemical reaction which enhances 
the reaction rates [33]. However, for the SFF design, the reac- 
tants enter the GDL by forced convection at the flow channel 
bends due to the corners as well as diffusion. In addition, there 
is also the under-rib convection [35,37,38] between adjacent 
flow channels. The forced convection and the under-rib convec- 
tion significantly enhance the local oxygen mass flow rates in 
the porous layers, so the flow channel width ratio has less effect 
in the SFF design. 
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Fig. 11. Effect of flow channel width ratio on cell performance. 
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Fig. 12. Oxygen mass flow rates at y = 11.5 mm on cathode GDL-CL interface 
at operating voltages of 0.3 and 0.7 V for triple loop flow fields with various 
flow channel width ratios. 


Figs. 12 and 13 show the oxygen mass flow rates and liq- 
uid water distributions at y= 11.5 mm on the cathode GDL-CL 
interface at the operating voltages of 0.3 and 0.7 V for the 
various flow channel width ratios. The liquid water distribu- 
tion is opposite to the oxygen mass flow rate distribution. At 
the higher operating voltage of 0.7 V, the oxygen mass flow 
rates and liquid water concentrations, are lower than those at 
the lower operating voltage of 0.3 V, with their values almost 
equal to each other for the three flow channel width ratios, 
so the cell performance is not dependent on the flow chan- 
nel width ratio. At the lower operating voltage of 0.3 V, as 
the flow channel width ratio increases, the oxygen mass flow 
rates increase slightly, while the liquid water concentrations 
decreases slightly, so the cell performance improves a lit- 
tle. 


4.4. Pressure drop losses 


The flow field design influences not only the cell performance 
but also the pressure drop in the fuel cell. Larger pressure drops 
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Fig. 13. Liquid water concentrations at y = 11.5 mm on cathode GDL-CL inter- 
face at operating voltages of 0.3 and 0.7 V for triple loop flow fields with various 
flow channel width ratios. 
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Table 4 

Estimated pressure drop losses at 0.3 V 

Flow field types AP (Pa) Ween (W m~?) Wp (W m~?) Wnet (W m~?) 
Single SFF with 7 bends 671 8245.52 14.80 8230.72 
Single SFF with 11 bends 1384 8820.99 30.52 8790.47 
Single SFF with 15 bends 1853 9027.16 40.88 8986.28 
Double SFF 927 7727.30 10.22 7717.08 
Triple SFF with n =0.4 1002 6519.46 7.37 6512.07 
Triple SFF with n =0.522 608 6564.96 4.47 6560.49 
Triple SFF with 7 =0.6 470 6577.37 3.46 6573.91 


in the fuel cell mean that more power is needed to pump the 
reactants. Thus, the pressure drop is a significant issue to be 
considered in choosing the flow field designs in addition to the 
I-Veey, curve. The cathode pressure drop has been converted to 
a power density using [21]: 


_ AP ‘A channel V 


(19) 
Atotal 


P 
where Wp represents the cathode pressure drop loss, AP the 
total cathode pressure drop in the fuel cell, Achannei the cathode 
cross-sectional flow inlet area, V the fuel velocity at the cathode 
inlet and Atotai represents the total reaction area. Table 4 lists 
the pressure drop losses for all the flow fields discussed in this 
paper. As the number of flow channel bends increases and the 
number of serpentine loops decreases, the total pressure drops 
and pressure drop losses increase, but the pressure drop losses 
are far less than the cell output power, Ween, for the miniature 
PEMEFCs due to their small flow channel cross-sectional areas 
and low reactant inlet velocities. Therefore, the pressure drop 
losses can be neglected. 


5. Conclusions 


Appropriate flow channel design is an effective means for 
improving cell performance. A three-dimensional numerical 
model was used to analyze the effects of various design param- 
eters for the bipolar plates on the cell performance and local 
transport phenomena for PEMFCs with the serpentine flow field 
design. The results indicate that: 


1. At high operating voltages, the various bipolar plate design 
parameters have little effect on the cell performance. 

2. At low operating voltages, the design parameters signifi- 
cantly affect the cell performance. For the single serpentine 
loop flow field, the cell performance improves as the num- 
ber of flow channel bends increases. The single serpentine 
flow loop field has better performance than the double and 
triple loop flow fields. The flow channel width ratio, n, has a 
relatively small effect on the cell performance, with the best 
performance for n = 0.6. 

3. Although the pressure drops differ for the various flow field 
designs, the pressure drop losses can be neglected for the 
miniature PEMFCs because the pressure drop losses are far 
less than the cell output power. 
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